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Abstract. The Fokker-Planck equation with diffusion coefficient quadratic in space variab- 
le, linear drift coefficient, and nonlocal nonlincarity term is considered in the framework of 
a model of analysis of asset returns at financial markets. For special cases of such a Fokker- 
Planck equation we describe a construction of exact solution of the Cauchy problem. In 
the general case, we construct the leading term of the Cauchy problem solution asymptotic 
C4_H , in a formal small parameter in semiclassical approximation following the complex WKB- 

JL_| Maslov method in the class of trajectory concentrated functions. 
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^J" : 1 Introduction 

00 

<^ • Studies of observations available at financial markets reveal empirically stated statistical regu- 

larities of asset price changes whose distributions are of broad interest providing more reasons to 
control of losses and gains in the course of quotations. Asset returns can be measured in terms 
of increments A£ = £(£ + At) — £(£) over the time interval At of a time series £(£) of quota- 
tions of an asset. The quantitative handling of high-frequency data of currency exchange rates 
(see pQ and references therein) have substantiated heavy tailed probability density functions of 
price increments A£ for small time steps At, asymptotic power law (or heavy tails) of return 
distributions different from the Gaussian distribution, volatility clustering, etc. 

Power law is considered to be well-established fact (see, e.g., [2]). Several semi-empirical 
mathematical models have been proposed to describe the non-Gaussian properties of asset return 
distributions. The most known of them are ARCH (auto regression conditional heteroskedastic 
models with time changing conditional variances) [3] and its modification GARCH |3], "trun- 
cated" Levy distributions [5], multiplicative noise models (see, e.g., [6]), multifractal cascade 
models with the so called HARCH effect [7|. 

The model proposed in pQ to describe asset return distributions uses the approach based on 
the Fokker-Planck equation (FPE) with variable coefficients. In some papers this equation is 
called the Fokker-Planck-Kolmogorov equation. The idea of the work pQ is based on similarity 
observation between asset return distributions and distributions of vortex velocity increments 
in turbulent cascades in hydrodynamics. This observation is verified by quantitative analysis 
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of distribution parameters that justifies an analogy between turbulent cascade in a stream of 
liquid and information cascade at stock markets. 

Following this semi-empirical analogy, we fix our attention on the power law of the probability 
density function (PDF) on one hand, and on the analogy with the turbulent cascade, on the 
other. For such a system, self-organization phenomena are the characteristics that could be 
mathematically expressed by an additional nonlinear mean drift term in the FPE [8]. 

In this work we modify the FPE with the diffusion coefficient quadratic in A£ considered 
in p] by introducing the nonlinear mean drift term. For such a FPE, we apply the approach 
developed in [9j [TOl [TT] which reduces the integrability problem for the nonlinear FPE to the 
problem for a linear equation associated with the original nonlinear FPE. We construct the exact 
solution of the nonlinear FPE whose coefficients have a special form, and for the coefficients of 
a general form we construct the leading term of solution asymptotic in a formal small parameter 
in semiclassical approximation in the class of trajectory concentrated functions. 

2 Formalism of the Fokker— Planck equation 

The formalism, proposed in pQ to describe distributions of stock market returns which asymp- 
totically obey power law, is based on the Fokker-Planck equation in the evolution variable At. 
The power law exponents are determined from the Kramers-Moyal coefficients founded in the 
framework of a phenomenological method developed in pQ. The FPE provides the information 
as to how the distributions on different time intervals At are correlated. In view of the analogy 
with the turbulent cascade the FPE is written in pQ in terms of logarithmic time-scale cascade 
proceeds. 

The FPE for the probability density function u(A£, r) of price increment A£ at a given 
logarithmic time scale r is derived in pQ in the standard way from the Chapman-Kolmogorov 
equation and is presented in the form 



3 T u(A£,r) 



- D V(AZ,t) + -^-,dW(AZ,' 



d(AO ^ s ' ' ' d(Atf 



u(A£,r). (1) 



Here d T = d/dr, and the Kramers-Moyal coefficients, the drift Z)W(A£, r) and diffusion 
D( 2 >(A{;,t) coefficients, are found directly from the financial data and the following approxi- 
mations were obtained 

£> (1) (A£,r) = -0.44A£, (2) 

D (2) (A£,r) = 0.003e" r/2 + 0.019(A£ + 0.04) 2 . (3) 

We assume the solutions u(A£, r) of equation ([1]) to be real smooth functions decreasing as 
|A£| — ► oo. Notice, that equation (pQ) can be written as a balance relation that makes obvious 
the conservation of the integral f_ u(A£,T)d(A£). Taking the initial function u(A£,0) to be 
normalized, we have 

X) 

u(Af,T)d(Af) = l, r^O. (4) 



The linear FPE (pQ), (|2|), ([3]) is solved numerically in pQ and quantitative agreement is obtained 
between the solutions and distributions constructed directly from the data of quotations, i.e. 
experimental data. 

A more detailed analysis of the power laws it(A£) ~ (A£)~( 1+At ' made in [12] shows that 
the exponents fi found in the framework of the FPE formalism are not in good agreement with 
experimentally obtained values. 
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Different ideas have been involved in studying of "heavy tails" of financial time-series dis- 
tributions. Mathematically, such nonlinear stochastic systems are described in terms of the 
Fokker-Planck equation with coefficients dependent on the distribution. We consider the feed- 
back caused by reciprocal effects of assets at the financial market. In accordance with [8] and by 
analogy between the exchange rates and the hydrodynamic turbulence [TJ [TJ] let us modify the 
FPE dl|) assuming that its coefficients D^ 1 ' and D^ 2 ' depend on the order parameter described 
by the first moment 

/oo 
A£u(A£,T)d(AO. (5) 

-oo 

In this work we consider analytical methods of solution construction for the following special 
case of the nonlinear Fokker-Planck equation (NFPE) 



d T u(Ati,T) 



(aAZ + xX u (T)) + -?—DM(AZ, 



<9(A£) v s " v " d(Af) v s ' 'd(Af) 



u(A£,r). (6) 



Here a and x are real positive parameters. The diffusion coefficient D^ 2 ' is taken in more 
general form as compared to (|3|), i.e. 

DM(A£, t) = e(/(r) + (aA£ + g{r)) 2 ). (7) 

Here /(t)(^ 0) and g(r) are real smooth functions, a is a real parameter, e is a real positive 
parameter. 

Market impact on price increment A£ of an asset with the probability density function 
u(A£,t) determined by ((6|), is modeled by the coefficient kX u {t). 

The motivation of the drift coefficient 

D^(A^T) = -aA^-xX u (r) (8) 

in can be taken from the stochastic methods [14] and empirical data handling [H [12] . In the 
Fokker-Planck equation 

d T u(A£, t) =d*t[- V A( (At t) + D^dAs] u(A£, r) (9) 

the function V(A£,r) can be treated as a potential of a regular external force Va£(A£,t) = 
<9a£^(A£,t) [14j . Let us assume that there are a number of mutually interacting assets on the 
financial market. Following the mean field concept [8], we modify the FPE ([9]) by substitution 

/oo 
W A dM-y,r)u(y,r)dy. (10) 

-oo 

Here the integral simulates the cumulative effect of other market assets on price increment A£ 
of a separate asset with a potential W(A£ — y,r) which depends on the difference A£ — y due 
to symmetry of mutual asset effects. 

In accordance with [H [12], where linear form of the drift Va^(A£, r) = aA£ gives close 
approximations of empirical data, we take modified drift coefficient (|10j) as — oA^ + x f (A^ — 
y)u(y,r)dy that results in -«A( — kX u {t), a = a — x, that, in turn, gives ([8|). 

To estimate the parameter x, one can take into account that xX u (r) in ([8]) is responsible 
for cumulative effect of all marketable assets on the price increment A£ of a given asset. It is 
assumed that in the stock market >cX u (t) can be estimated with the help of an appropriate 
market index, as it characterizes the market as a whole. In the exchange operations (FOREX 
market), where indices are not used, such a characteristics could be deduced from the correlation 
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matrix of currency exchange rates. Anyway, the estimation of the parameter x is a subject of 
a special study. 

The integrability problem for the FPE in applications is practically solved numerically. As 
a rule, analytical approaches are based on special classes of functions in which solutions are 
constructed. A more detailed review of solution methods for the NFPE ([6]) can be found in |llj . 

Using the approach considered in [IT] we obtain the evolution operator for equation ([6]) with 
a constant diffusion coefficient D^ 2 ' = const. This result is used to construct exact solutions 
of ([6]) for the functions f(r) and g(r) of a special form. Following [10] for these functions of 
general form, we construct asymptotic solutions of ([6]) in semiclassical approximation in a small 
parameter e — > in the class of trajectory concentrated functions. 

3 Evolution operator for the FPE with constant diffusion 

Let us assume that in the NFPE ([6]) the diffusion is constant, D^ 2 '(Ax,t) = e. Setting A£ = x 
for convenience, we rewrite © as 

d T u(x,r)= [d x {ax + xX u {t)) + ed xx ]u(x,r), (11) 

where d x = d/dx, d xx = d 2 /(dx) 2 . Equation (jlip plays the basic role in our study of the 
nonlinear Fokker-Planck equation ([6]) and in the linear case (x = 0) it describes the Ornstein- 
Uhlenbeck process [IB] (see also [H]). 

Let us take a fixed value r = s as the initial time and consider the Cauchy problem 

u(x, s) = ^(x) (12) 

for equation (jlip in a class of functions of the Schwartz space S in the variable x for r ^ s. In 
accordance with §%§ we can set 

/+oo 
^(x)dx = 1 
-oo 

without the loss of generality then 

u(x,r)dx = l, t ^ s. (13) 

For the first moment ([5]) of the function u(x, r), 

/+oo 
xu(x,r)dx, (14) 

-oo 

from (jlip and (jT3j) we immediately obtain 

X u {T) = -(a + x)X u (T). (15) 

Here X u (r) = dX u (r)/dT, then 

X u (t) = X 7 exp [-(a + x)(r - s)] , (16) 

/+oo 
xj(x)dx. (17) 

-oo 

Therefore, the first moment X u (t) is determined by the initial distribution 7(2; ) and does not 
require the solution of equation ([fl]) . The substitution of (fTBP into (fTT|) actually linearizes the 
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NFPE. Such a property of the NFPE with coefficients depending on the first moment X u (t) is 
related to the fact that this equation is close to the linear one (see |11] for more details). 

Let us present the solution of the Cauchy problem (jlip and (|12p as action of the evolution 
operator U(t,s, •) on the initial function 7(2;). Following [11] we obtain 

/+00 
Giin(r, s,x - X u (r),y - X J )j(y)dy. (18) 

-00 

Here Gn n (r, s, x, y) is the Green function of the linear FPE resulting from (jllj) with the assump- 
tion x = 0, 

a(r-s) 



( x - e -(r-s) y y 



4esinha(r — s) 



ef0-«) 
Gn n {T, s, x, y) = - exp 

^ sinh a{T-s) 

Notice that the solution of the Cauchy problem (jllj) and (|12p can be immediately generalized 
for the multidimensional FPE. 



4 Nonlinear FPE with quadratic diffusion 

Consider the NFPE 

d T u(x,r)= [d x (ax + x(3(T))+ed x (f(T) + (ax + g(T)) 2 )d x ]u(x,T), (19) 

where coefficients /3(r), /(r), and g(r) in general case can depend on X u (r). 

With the obvious notation A£ = x, equation (fl9j) takes the form ([6|), ([7]) when /3(r) = X u (t), 
and /(r) and g(r) do not depend on X u (t), and, therefore, (fT9j) generalizes (|6j), (J7]). 

From (Tl3j) . (fTlj) . and (fT9l) by analogy with (fT5j) we immediately obtain 

^«(r) = -(a - 2a 2 e)X u {r) - x[3(t) + 2aeg{r). (20) 

Similarly to equation (jlip . the Cauchy problem (j!2p for the NFPE (j!9p induces the Cauchy 
problem (j20j) and 



X u (s) = X 7 (21) 

determining X u (r) independently of the solution of the FPE (j!9p , where X 7 has the form (|17|) . 
When p(r) does not depend on X u (t) and /?(t) = X u (t), we have the following explicit expres- 
sion for X u (t) from ^ and (|2Ijl 



J s 



The solution of equation (j20p can be found also for a different form of dependence of g(r) 
on X u (t). In any case, finding of the moment X u (t) separately from solution of the nonlinear 
FPE and substitution X u (t) into equation (|19p actually linearizes this equation with respect to 
thePDFu(x,r). 

In order to derive exact and approximate solutions it is helpful to introduce new variab- 
les (y,r), 



V 

a 



i log (ax + g(r) + V/(r) + (ax + 5 (r)) 2 ) , (22) 



G 
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such that in the new variables the diffusion term becomes a constant. In the new variables 
equation (|19|) takes the form 

q(y,r) 



d T u{y,r) 



a 



-d v + ed v 



Here the following notations are used: 



u(y,r) 



xftr) - - 9 (r) - 9 ^1)- M e - 
a a J a 



'HI 



O 



+ (- + ae) A { _ } {y,r), 



A {±) (y,r)=e a y±f(r)e 



<m 



(23) 

(24) 
(25) 



Let us impose special constraints on the functions /3(r), g(r), and /(r) to reduce ([23]) to 
a solvable form and to construct exact solutions. 



4.1 Exact solution 

An exact solution of equation (|19p can be found if the functions /(r) and g(r) obey the conditions 



/(r) + 2oc/(t) = 0, g(r) + a 5 (r) - ax(3{r) = 0, 



a 
c = — h ae. 
a 



(26) 



In this case f(r) = f(s)e 2ac(j ' s \ g(r) = g(s)e alyT s '+a>c/ T e a ( T ^/3(£)d£, and equation 
takes the form 

[— d T + cd y + a + edyy] u(y, r) = 0. 



(27) 



By changing the variables (y, r) — > (z, r), where 2 = y + cr and taking v(z, r) = u(y, r)e QT , we 
reduce (|27|) to the standard diffusion equation 

(28) 
(29) 



[-a T + e5 22 ]u(2;,r) = 0, 
whose Green function is well known 

1 (z-z') 2 



G st {T,S,Z^) 



_g 4e(r — s) 



•y/47re(T — s 

where a fixed value r = s is taken as the initial time. 
Then, for equation (|27p we have 

= a (T- S ) ( 1 



G{T,s,y,y') 



and for equation (|19 



^4^6 (t — s) 



exp 



G(t, s,x,x) 



a(r-s) 



aA 



exp 



7re r — s 



4e(r - s) 



1 

4e(r - s) 



[y - y' + c(t - s)]' 



1 ax + 5 (r) + V?(r) + (ax + 5 (r)) 2 
log 



+ c(r - s) 



a ax' + g(s) + ^//(s) + (ax' + g{s)) 2 

(30) 



when conditions (|26p are satisfied. Notice that the Green function ()30|) can be considered as 
kernel of the correspondent evolution operator of the NFPE (|19p . i.e. 

/+oo 
G(t, s, x, x')(/(r) + (ax' + ^r^r^xOdx'. 
-00 

The evolution operator U is nonlinear because the functions /3(r), /(r) and g(r) entering the 
Green function ([30]) depend on the function u(x' , s) = j(x') with the operator U acting on it. 
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4.2 Semiclassical approximation 

For equation (f23|) with the coefficients /3(r), /(r), and g(r) that are generally dependent on r 
but independent of e we construct a leading term of formal solution of the Cauchy problem 
and the evolution operator asymptotic in small parameter e (e — > 0) in the class of trajectory- 
concentrated functions. 

Let us introduce class V^(Y(t, e)) of functions 

where Ay = y — Y(t, e), y is given by (|22|) . and ^(r, e) is a functional parameter of the class, 
regularly depends on e and is to be determined. The function </?(£, r, e) belongs to the Schwartz 
space, smoothly depends on £ and r, decreases with £ — > oo, and regularly depends on e for 
e — > 0, i.e. it can be expanded in the Taylor series around e = 0. 

Functions of the class V^(Y(r, e)) are normalized with respect to the norm 



JR 1 



By analogy with |1Q|. I15j let us call solutions of the class VJ(Y(t, e)) the trajectory concentrated 
solutions. Restrict ourself with the case when Y(t, e) does not depend on e, Y(t, e)= Y(t). We 
are interested here in solutions of equation (j23[) localized in y in a neighborhood of the curve 
y = Y{t) for each fixed e G [0, 1[ and rel 1 . 

For the functions belonging to VJ(Y(t)), the following asymptotic estimates can be imme- 
diately obtained 

||(Ay) fc *||(r, e ) _ k/2 \\(ed y ) k n(r,e) _ k/2 

\\H(r,e) U[C h ||*||(r,e) " U{ * h 

\\e(d T + Y(T,e)d y )<S>\\(T,e) 



m\(r,e) 
y-y{r) 



0(e). 



Notice that £ = - — -X-L has the estimate £ ~ 0(1). 

Transformation from (y, r) to the variables (£, r) reduces equation ([23]) to the form 



Ltp(£,T,e) 



a l j. (v( \ . g(v^i+z(i)ii) ^ x fl L P 



p(£,T,e)=0, (31) 



where notations Q24|) . (j25|) are used. 

Let us expand coefficients of equation (f3T|) and the function </?(£, r, e) into power series in y/e, 
e.g. 

<^,r,e) = <p<®(Z,T) + vV^CM + • • • . (32) 

Equating coefficients at equal powers of (y/e) < 2 , A; = 0, 1, 2, . . . , in (I27p we obtain the following 
equations determining asymptotic solution of equation (|3ip 



r(r )+ ;°'yff-;). = o, (33) 

A (+) (F(r),r) 

i(0)^ 0) (e,r)=0, (34) 

L {0) <pW(a,T) = -L {l) ^ \!i,T). (35) 
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Here we use the notations 



L(o) = -d T + % + f -R(t)<% + a, 
L {1) = (Qo(r)e + Qi(r)) d^, 



R{t) = 
Qo(r) 



A i+) {Y(r),r) 



2A {+) (Y(t),t) L 
aA(_)(y(r),r) 

" %(nr),r) 



f(r)e- aY ^ + aA {+) (Y (t), r) + aF(r)A ( _) (y(r), r) 
/( r)e -^M + aA( _) (Y(r), r) + aY{r)A {+) {Y(t),t) 
aA H (Y(r),r) 



R(t), Q 1 (t) 



A (+) (Y(t),t) 



q (y,r) = q(y,T,e)\ e=0 = 2 fx/3(r) - -g(r) - ^-) - f M e -*v + ^A^t). 

\ a a J a a y ' 

It is obvious that the construction of the asymptotic solution is reduced to the solution of 
equation ([M|) . 

By changing dependent and independent variables in ([34"|) according to the formulas 



(0) 



ti,T) = e a ^v(z,r f ), t'= / h 2 ( V )dri, 



V 



z = h(r)£, h(r) = h(s) exp / R(rj)dq , 



(36) 



we reduce equation ([3"4"|) to equation ([2S]) with e = 1 whose Green function is given by (|2"9"j) . i.e. 



G st (r',0,z,z') 



V / 47TT 



■.e 4t' 



(37) 



Notice that in (136j) r = s corresponds to r' = 0. Therefore, the solution of the Cauchy problem 
( — d r i + d zz )v(z,r / ) = 0, v(z, 0) = i>o(z) given by the Green function ([37]) reads 



(38) 



u(z,t') = / G at (r / ,0,z,« / )uo(« / )d^ 
and for the problem ( — d T > + d zz ^u(z, r') = F(z, r'), u(z, 0) = 0, we have 



u(z, r') = d V G st {T' - V , 0, *, z')F{z', rj)dz' . 

JO J -co 

Let us now find the asymptotic solution of the Cauchy problem 

<p(£,s) = 7(0 

for equation (I3ip in the form (|32|) . To this end we change variables from (z,t') to (£,t) in (J38l 
according to (|36|) and obtain the solution of the Cauchy problem ip^°'(^,s) = 7(£) for equa- 
tion (El) 



(39) 



(40) 



</»(£, r ) = / e «(— )G st (r'(T)Ah(r)C,h(s^'MaK^' 



-OO 

oo 



= a(r-s) 



4:7: J T h 2 (r])dr] 



cxp 



(MT)g-ft(*)O s 

*JJh 2 (r,)dr, 



y(?)h(8)d£. 



(41) 
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Here the relations z' = h(s)£' , j(£') = vq(z') = vo(h(s)£') are used. Similarly, in accordance 
with (|39|) for the function </? (£> T ) determined by ([35]) we have 

/t /*oo 

drj / e a ^G st (T'(T)-T'( V ),0,h(T)Z,h( V )aF(Z',v)Hv)dZ', (42) 

J — OO 

where F(h( V )e ,r'(rj)) = e-^-^h~ 2 (r,)F(e,v), HM = -L (1) ^°)(£,t), 

{h( T )z - h{ v )Z') 



e a(T-r,) 

<P ( '(£i T ) = d V i ^=exp 



4j;h*(y)dy 



F(?,r,)h( V )dt'. (43) 



4nrfih?(y)dy 

Thus, the function ([32]) determined by expressions ([33]) . ([3T]) . ([32]) gives the asymptotic 
solution of the Cauchy problem (I3ip . (|40j) accurate to 0(e) in the following meaning. Denote 
the action of the operator L defined by ([3"T]) on the function <£ (£ ! ' r ) + V^ (£j t ) as 

S (2) (£,T,e) :=L(^(e,r) + ^ (1) (e,r)), (44) 

then 

Il5 (2) ll(r,6) 



lvll(r,c) 



O(e). (45) 



Similarly to (|18|) . we can consider the solution of the Cauchy problem as action of the evolution 
operator of the operator L on the function 7(£). 



5 Concluding remarks 

Asset interactions at the market are naturally supposed to be complex and nonlinear. The ap- 
proach discussed in the work is intended for the investigation of the asset market by using simple 
nonlinear model directly generalizing known models based on the linear FPE ([I]). The asymp- 
totics obtained, equations (J32J) . (J4TJ) . (f43l) . can be considered as a necessary step in studying of 
nonlinear properties of the asset market. 

The approach provides the way to analyze effects of nonlinearity on the distribution dynamics 
of asset price increments at the stock market or currency exchange. Such an analysis necessarily 
implies the statistical processing of corresponding empirical data to estimate values of the model 
parameters, which is beyond the scope of the present work and is a subject of a special research. 

Following [H [12j [13] , we use quadratic approximation for the diffusion coefficient and li- 
near drift ([8]). Our formalism can be also applied to the NFPE with more general dependence of 
coefficients on A£. The quadratic approximation of the diffusion coefficient ([7]) is a generaliza- 
tion of the constant diffusion that results in non-Gaussian probability density functions. These 
PDFs are asymptotic solutions of the NFPE (|19|) in semiclassical approximation in explicit form 
given by (|32|) . (J4T1) . (|43l) accurate to O(e) in accordance with ([31]) . ([35]) . The construction of 
semiclassical solutions is based on exact solution of the NFPE with constant diffusion (|11[) . 

Preliminary calculations show that for the initial function 7(3;) having the power law asymp- 
totics ("fat tails") the distribution keeps the non-Gaussian shape during its evolution. 

In conclusion we would notice that since the stochastic properties of asset price increments are 
studied effectively by numerical simulation, it is of interest to formulate stochastic differential 
equations with the nonlinear feedback related to the NFPE considered above. 
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